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1. Introduction 
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Cluster formation in strongly interacting Fermi-systems is one of the most interest- 
ing subjects in many body physics. So far, in condensed matter, only the study of 
formation of two body clusters has been studied in a wide perspective in the con- 
text of pairing, that is superfluidity and superconductivity (see e.gP in the context 
of cold atoms). However, in nuclear physics, due to the existence of four different 
fermions, all attracting one another, the scenario of cluster formation is much richer. 
Of course, there also exist Cooper pairs of neutrons and protons, and to a lesser 
degree proton-neutron pairs, but in addition there exist tritons, helions (trions) and 
a-particles (quartets), to cite only the lightest ones. For instance, the formation and 
condensation of a-particles is presently in the focus of our studies. However, adding 
to a gas of a-particles two neutrons per a-particle, at some higher excitation energy, 
this system may be transformed into a gas of tritons. For example, one may imagine 
that 12 Be is composed out of four tritons at sufficiently high energies. Also nucleons 
themselves are trions when seen as bound states formed out of three quarks. Then 

"This work is part of an ongoing collaboration with Y. Funaki, H. Horiuchi, A. Tohsaki, T. Yamada. 
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the study of the transition from a Fermi gas of quarks to a new Fermi gas of trions, 
i.e. nucleons is interesting. 

Recently in cold atoms physics the trapping of three different fermions has been 
achieved by two experimental groups and the formation of a gas of trions seems 
possible ®. May be in the future, one will even capture four different fermionic 
atoms and study quartet condensation. Theoretical work on trions and quartets 
already has appeared in quite a number (see articles in S and references therein) 
and it seems that the field will rapidly expand in the future. The theoretical works 
mostly have been performed in one dimensional (ID) models. In this work, we will 
try to develop a many body approach for the treatment of a gas of clusters of 
fermions valid in 3D. 



2. Alpha-particle Condensation 

The condensation of a-particles can be treated in analogy to the pairing case. To 
this end one has to consider the equation for the four fermion order parameter ^1234 
(summation convention is used) SI 

1 — m — 7J2 

V1234 = — ; fi2i'2' v I / i'2'34 + permutations (1) 

V - £1234 

where [i is the chemical potential and £1234 = £1 + £2 + £3 + £4 where the q's are the 
single particle energies. Besides the matrix elements of the interaction um'2', the 
most important ingredients are the single particle occupation numbers n% — (cfci). 
They, in principle, have to be calculated from a Dyson equation with a mass operator 
containing the quartet condensate. Equation ([T]) then constitutes the quartet "gap- 
equation" si. 

A lowest order expression for the single particle mass operator entering the 
Dyson equation can be obtained in the following way where the four particle T- 
matrix T4, in the one pole approximation corresponding to (TTJ), is convoluted with 
three uncorrelated hole lines. 

£(1,V)= T 4 (1234, 1234; n 4 )G3(234,fi 4 - 2*0 (2) 

234, fi 4 

where we used summation over the Matsubara frequency and the uncorrelated three 
body Matsubara Green's function is given by 



G°(234, z v ) = ( 1 --ft)(l-/ 3 )(l-/4) + /2/3/4 (3) 
z v — £2 — £3 — £4 

However, in order to simplify in a first step, we linearize and replace, at fi- 
nite temperature, the m by the Fermi-Dirac occupation numbers Uj — > f% = 
[1 + e( £i ~ ,1 '' T ]~ 1 . Still (H|) is very difficult to solve, since it is a four body problem 
rendered more complicated by the presence of the Pauli blocking factors. However, 
for clusters involving more than two fermions, and certainly for strongly bound 
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quartets like the a-particle a projected product ansatz^lfor ^1234 is known to be 
a good approximation: 

* 1234 -> S(K - ki - k 2 - k 3 - k 4 ) ( /) O (k 1 ) ( /) O (k 2 )0 o (k 3 )0 o (k 4 )x ST (4) 

where K is the total c.o.m. momentum and <fro(k) is a mean field OS 1 - wave function. 
Spin and isospin are taken care of by x ST with S = T = for the a-particle. We will 
not consider it further. Inserting Q into ([1]) (with K = 0, for condensation), one 
straightforwardly obtains a Hartree-Fock (HF)-type of equation for 4>q which can 
be solved for various critical temperatures T = T c , while \i is determined from the 
particle number condition N/V = 4/ ^h) 3 f( £ p)> see elsewhere ® for the detailed 
calculation. For simplicity we only take a spin-isospin averaged separable interaction 
Vl2V2 , = Xw(^^-)w(^^-)S{k 1 + k 2 - ki - k 2 ) with the form factor w{k) = 
e - fe2 / fe o p Xhe two open parameters A = -991 MeV fm 3 and fco = 1.43 fm -1 are 
adjusted that energy (-28 MeV) and radius (1.71 fm) of a-particle come right. In 
Fig. 1 we show the result for T c as a function of p, and the density n where the free 
gas Equation of State (EOS) has been used to relate n and fi. It is very remarkable 
that the obtained results for T" well agree with a direct solution of (1) using a 
realistic NN forced These results for T C Q are by about 25 percent higher than the 
ones of our earlier publication SI. We, however, checked that the underlying radius 
of the a-particle in that work is with 2 fm larger than the experimental value and 
that T" decreases with increasing radius of a. Furthermore a different variational 
wave function was used in SI. 

In Fig. 1 we also show the critical temperature for deuteron condensation. In 
this case we take A = -1305 MeV fm 3 and fco = 1.46 fm -1 to get experimental 
energy (-2.2 MeV) and radius (1.95 fm) of the deuteron. It is seen that at higher 
densities deuteron condensation wins over the one of a-particles. This stems from 
the fact that Fermi-Dirac distributions in the four body case, see ([1]), can never 
become step-like, as in the two body case at zero temperature, since the pairs in an 
a-particle are always in motion. As a consequence, a-condensation only exists as a 
BEC phase and the weak coupling regime is absent. 

The mean field ansatz of course, simplifies the problem of quartet condensa- 
tion enormously. In the future it, therefore, should be possible to solve the quartet 
gap-equation |(5J) with f4} fully self- consistently. At zero temperature this will corre- 
spond (approximately) to the minimization of the energy with a four body coherent 
state 

|tf Q > ~ e* £"34 *(1234)c+e+c+c+ (5) 

The relation between ^1234 and <?!>(1234) is not trivial but can be established in 
an analogous way to BCS theory. 

Condensation phenomena can have precursor signs in finite systems, even small. 
We know that, e.g., from nuclear pairing. Evidence for a-particle condensation exists 
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in light self conjugate nuclei. For example the so-called Hoyle state at 7.65 MeV in 
12 C turns out to be a product state of three a-particles in identical OS" orbit to over 
70 percent Indications from theoretical results have been put forward recently 
that the 6-th + state at 15.1 MeV in 16 is a four a-particle condensate stated If 
this prediction is further verified experimentally and theoretically, then a-particle 
condensation very likely is a general phenomenon in nuclear systems Of course, 
for the description of the condensation of a very small number of a- particles a 
number conserving variant of ([5]) must be employed. Also the c.o.m. motion is not 
described by plane waves as in (TjQ) for homogeneous systems but the c.o.m. orbitals 
obey themselves a mean field equation. The corresponding ansatz, therefore, is 



n,r 2 , . 



.4[0(1234)< ? (>(5678). 



(6) 



1234 



with A the antisymmetriser and 0(1234) = <I>o(Ri234)v3o( r ij) where R 
(ri+r2+r3+r4)/4 is the c.o.m. coordinate and = Vi—Yj stands for all the possible 
combinations of relative coordinates of the quartet. Of course, at the end the depen- 
dence on the total c.o.m. coordinate of the system Xq = [R1234 + R5678 + ■■■]/n a , 
with n a the number of quartets has to be eliminated , that is one has to project on 
good total c.o.m. momentum zero. In the case of a Gaussian ansatz for $o(R) and 
ipo(rij), this is easy to write down and we obtain 



lintrinsic] 
y J 



(7) 



with 



V 



-(R?234)7BV(R5 G 678) 2 /B 2 



(8) 



_ s 2 /h 2 -s 2 /b 2 



(9) 



where Rf 2 34 = R1234 - X G , and s 2 234 = (n - r 2 ) 2 + (ri - r 3 ) 2 



alpha (rms radius: 1.71fm) 
deuteron (rms radius:!. 95 fm) 
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Fig. 1. Transition temperature of alpha and deuteron as function of chemical potential (left) and 
density (right). 
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We see that the a-condensate wave function consists in the anti summarized 
product of two interdependent mean field parts, one concerning the intrinsic part of 
each a-particle and one concerning the c.o.m. motion of the a-particles. Each part 
is a mean field wave function with the c.o.m. part eliminated. The wave function 
(|6|) is, therefore, totally translational invariant. This type of wave function is known 
in the literature as THSR wave function^. It depends on two parameters B and 
b which in ^ were taken as Hill- Wheeler coordinates to quantize the corresponding 
energy surface Eo(B, b). For B 3> b one can neglect the antisymmetrizer in ([6|) and, 
as we see from ((5]), the wave function becomes a pure product state of a- particles, 
i.e. a condensate. On the other hand for B = b ([6]) is equivalent to a pure Slater 
determinant The ansatz (j6])-([8]) gives excellent results for all observables of the 
Hoyle state with the Pauli principle only acting mildly^. It also predicts similar a- 
condensate states for 16 0, 20 Ne, ... close to the corresponding threshold energies for 
a-particle break up ®. For the case where in nuclei many a-particles existed (e.g. 
40 Ca with 10 a-particles), one also must employ the coherent state formulation, 
since the explicit antisymmetrization becomes very difficult. 

3. Trions 

The case of a gas of trions has not been considered very much in nuclear physics. 6 He 
has an excited state at ~ 12 MeV which can be interpreted as mainly consisting 
out of two tritons One can speculate that, e.g. in 12 Be there exists a state 
around 30 MeV excitation energy with four tritons. Since tritons are fermions, the 
corresponding wave functions of the c.o.m. motion will develope nodes as in the 
usual shell model with OS, OP, etc. orbitals. The THSR wave function for a few 
tritons may have the form 

\nt) ~ ^(123)^(456)....] (10) 

with 

0(123) = e- R »/ B Vt(r I3 ) (11) 

with definitions analogous to the ones of the quartet case. As already mentioned, 
the antisymmetrizer will make out of the product of c.o.m. Gaussians a Slater de- 
terminant, whereas the intrinsic wave functions of the tritons may stay essentially 
unaltered. 

In the nuclear case, a gas of tritons will contain only a very small number of 
such clusters. Considering nucleons as trions formed by three quarks, the number of 
trions can be considerable, that is the number of nucleons in a nucleus or in a Heavy 
Ion Collision. In the case of cold atoms where, as mentioned, recently trapping of 
three different fermions has been achieved the number of trions may go in the 
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ten thousand. In such cases the trion wave function \nt) of above cannot be handled 
directly and we must elaborate some approximate scheme. To put a trion creation 
operator in the exponent like for the quartet is not possible, unless one introduces 
Grassman algebra, since we deal with fermionic clusters. One way to proceed is via 
the well known Equation of Motion (EOM) method , in analogy to Self Consistent 
Random Phase Approximation for two particles In this case one first defines in 
the present case a three particle creation operator 



Q^2>i23^ + 23 ( 12 ) 
123 

with C123 = c i~ c 2 c 3 /[(I — Tl i)(l — n 2)(l — ^3) + nin 2 n 3 ] 1 / 2 the product of three 
fermionic single particle creators. The equation for the amplitudes x can be obtained 
by minimizing the following energy weighted sum rule: 

_ ([Q u [H,Q+]]+) 

1 (IQuQth) (13) 

This leads to a non linear hermitian secular equation for the x's. The nonlinear- 
ity stems from the fact that in the double commutator in ||13|) appear three-body 
correlation functions of the type (c + c + c + ccc). The relation between Qf and C± 23 
can be inverted and the three body correlation function be expressed by the x's 
with the help of the relation Qt\0) — 0. However, two body and one body cor- 
relation functions appear as well. Though approximate formulas expressing those 
lower rank correlation functions by the x-amplitudes may be derived, thus obtaining 
a fully self consistent system of equations, the procedure seems extremely heavy. 
A strong simplification could consist in factorizing all correlation functions into an 
antisymmctrized product of single particle correlation functions (single particle den- 
sity matrices) and express them via x-amplitudes in a similar way as we discussed 
above for the alpha particle condensat ion . This procedure is generally known as the 
so-called renormalized RPA, see, e.g., UJJ 

On the other hand, if in the gas the trions stay more or less compact, one may 
try a similar mean field ansatz as for the a-particles, inspite of the fact that a 
mean field description of only three fermions may not be as good as the one for 
a quartet. However, in the case of three colors, the three fermions can occupy the 
lowest OS" level of the mean field and with a projection on good total momentum 
the description may still be quite reasonable. For instance, this may be a quite 
attractive approach for the constituent quark model where the three quarks in a 
nucleon are bound by a harmonic oscillator field. Let us, therefore, make a THSR 
ansatz also for the trion wave function for the case of homogeneous infinite matter 



X123 -> e KRl23 5(K - ki -k 2 - k 3 )^(k 1 )^(k 2 )^(k 3 )Ci23 (14) 
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where (123 takes care of spin, isospin, and color indices. With (TT4")) we can insert 
this trial wave function into (|13[) and vary with respect to y(k). Since the trions are 
fermions, we also must take care of the fact that the trions should carry different 
c.o.m. momenta K and, therefore, in (|12p a summation over the c.o.m momenta 
must be contained. The case of trion formation has been studied theoretically in 
solving exactly ID-model Hamiltonians with three colors 5. It would be interesting 
to compare the results of above approach to these exact solutions. 

In the case of nucleons, one pair of quarks is forming a strongly correlated di- 
quark state. It seems reasonable to approximate the di-quark by a boson. Then 
the transition of a gas of quarks and di-quarks can be considered as a mixture of 
fermions and bosons. The nucleons then constitute a bound state between a fermion 
and a boson. Interesting physics has recently been found in treating the fermion- 
boson scattering in the background of a gas of bosons and fermions (the one pair 
approximation). Most interestingly the fermion-boson pair seems to form a stable 
pair even for infinitesimal attraction, quite similar to what happens for the formation 
of Cooper pairs in a two component Fermi gas. However, the boson- fermion pai rs 
have, of course, Fermi statistics, thus forming a new Fermi gas of composites 

4. Pairs 

In the pairing case, the focus is on the BEC-BCS transition. Researches of BEC- 
BCS crossover have been done not only for the condensed matter 03 but also for 
the nuclear matter M The application of the Nozieres Schmitt-Rink approach^ 
is fairly standard by now, taking care of the purely bosonic (strong coupling) limit 
of the pairs and their c.o.m. motion in the Bose-Einstein distribution. However, 
in the strong coupling limit, the pairs behave as (composite) particles and create 

I -j H I 

their own mean field. This effect is not taken into account in the work of as 
discussed by the authors themselves. In order to include this effect, we can derive 
two coupled mean field equations, one fermionic, the other bosonic (fermion pairs). 
This is achieved by the minimization of the following generalized sum rules 

e k = ([q k [H,q+]] + )/([q k ,q+] + ) (15) 

with q^ = u k c£ — ^fcCfc being the standard Bogoliubov transformation of fermions. 
The corresponding bosonic equation is obtained from 

Ek, v = ([Qk,„, [H, <&,„]]>/<[<&£,*, <&,„]) (16) 

where 

Qk,u = E - i&^J/a - n kl - « k2 ) 1/2 (17) 

kik 2 

represents a Bogoliubov transformation for fermion pairs. 
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A momentum conserving delta function between K and ki , k2 is implicit and spin 
and isospin indices have been suppressed. The expectation values are with the cor- 
related ground state defined (approximately) by Q\0) = 0. 

Equations (|15[) -(|lT p can be generalized to finite temperature in introducing cor- 
responding Gorkov equations. The fully selfconsistent set of equations for ampli- 
tudes u, v, X, Y can be worked out (see and references therein) and their solu- 
tions give very promising results in non-trivial model cases HSJ. They constitute a 
fully selfconsistent extension of the Nozieres Schmitt-Rink approach 

5. Conclusions 

Selfconsistent many body approaches for gas phases of quartets, trions, and pairs 
in multicomponent Fermi-systems have been established. They are applicable to 
realistic 3D systems accounting for strong cluster phenomena. 
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